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qh| machine, and continues to be of enormous benefit in the classical the- 

ory of computation. We propose that quantum computation, like its 
classical counterpart, may benefit from a version of the lambda calcu- 
lus suitable for expressing and reasoning about quantum algorithms. 
' In this paper we develop a quantum lambda calculus as an alternative 

model of quantum computation, which combines some of the benefits 
of both the quantum Turing machine and the quantum circuit mod- 
els. The calculus turns out to be closely related to the linear lambda 
calculi used in the study of Linear Logic. We set up a computational 
model and an equational proof system for this calculus, and we argue 
that it is equivalent to the quantum Turing machine. 



AMS Subject classifications: 81P68, 68N18, 68Q10, 03B70 
Keywords: Quantum Computation, Lambda Calculus, Linear Logic, 
Models of Computation 
Brown preprint: BROWN-HET-1366 



1 



1 Introduction 



Currently there exist two main approaches to the theory of quantum com- 
putation: The quantum Turing machine, introduced by Benioff and Deutsch 
[HE], and the quantum circuit model, introduced by Deutsch [Hj. These two 
approaches were shown to be essentially equivalent by Yao |1] . 

The quantum Turing machine provides a fundamental model of quantum 
computation that may be regarded as a baseline for defining universality. 
However, reasoning about Turing machines can be a cumbersome process, 
requiring word-at-a-time thinking while keeping track of complicated ma- 
chine and tape states. Turing machine programs do not satisfy a simple 
algebra. 

For this reason, the quantum circuit model is more popular in the practi- 
cal investigation of quantum algorithms. Quantum circuits are visual, com- 
positional, and may be manipulated algebraically. However, no single finite 
quantum circuit is universal. Indeed, Yao's proof of Turing equivalence relies 
on the concept of uniform circuit families generated by classical computation 
jH |3] . To define what we mean by such a circuit family, we need to rely on a 
separate model of classical computation not described by any finite quantum 
circuit. 

In classical computation, the lambda calculus provides an alternative 
computational model, equivalent to the Turing machine, which continues to 
be of enormous utility in the theory of computation, in mathematical logic, 
and in the study of computer languages and their semantics jHl 13 El 13 UD] • 
Due to its simplicity and expressive power, the lambda calculus has been 
used as the basis of several powerful computer languages, including Lisp, 
Scheme, ML and Haskell dH |I2l 1131 IH • 

In this article, we propose that quantum computing, like its classical 
counterpart, may benefit from an alternative computational model based on 
a version of the lambda calculus suitable for expressing and reasoning about 
quantum algorithms. We develop such a calculus, which turns out to be 
closely related to the linear lambda calculi used in the study of Linear Logic. 
We set up its computational model and equational proof system, and argue 
that the computational model is equivalent to the quantum Turing machine. 

The quantum lambda calculus combines some of the benefits of both 
quantum circuits and the quantum Turing machine. The quantum lambda 
calculus describes functions that may be composed and manipulated alge- 
braically, like quantum circuits. Programs can be algebraically transformed 
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into equivalent programs, and one can solve equations whose unknowns are 
programs, in much the same way as one solves equations in high school al- 
gebra [inj. Unlike quantum circuits, the quantum lambda calculus provides 
a unified framework that is universal for quantum computation without the 
need to rely on a separate model of classical computation. 

In a practical vein, we show how various known quantum algorithms 
may be expressed as simple programs in the lambda calculus. Indeed, the 
calculi described in this paper may be used as a programming language for 
prototyping quantum algorithms. In fact, the algorithms exhibited in this 
article were transcribed into Scheme for testing. The simulator, which was 
also written in Scheme, is available upon request from the author. 

Since the first version of this paper was written, some progress has been 
made by the author in devising a typed version, with accompanying deno- 
tational semantics, of a fragment of the quantum calculus described here 

2 The classical lambda calculus 

We begin by providing a reasonably self-contained introduction to concepts 
and constructions in the classical lambda calculus that will be used in the 
rest of the paper. The intended audience for this section includes physicists 
and general computer scientists. The expert may skip this section and refer 
back as needed. 

The classical lambda calculus may be regarded both as a programming 
language and as a formal algebraic system for reasoning about computation. 
It was originally introduced by Church in the study of the foundations of 
mathematics f7| ^] . Church postulated that it provides a universal model 
of computation, which was later shown by Turing to be equivalent to the 
Turing machine |19j . 

As a formal system, the lambda calculus has axioms and rules of inference, 
and lends itself to analysis using the language and tools of mathematical logic. 
Computation may be regarded as guided deduction in this formal system. 
This provides a directed form of equational reasoning that corresponds to 
symbolic evaluation of programs via a sequence of algebraic simplifications 
called reductions El EDI • 

The syntax of the classical untyped lambda calculus A is as follows: Ex- 
pressions (also called terms) are constructed recursively from variables x, y, 
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. . . , parentheses, spaces, the period, and the symbol A, according to the 
grammar of figure d 



t ::= 


terms: 


X 


variable 


{\x. t) 


abstraction 


{t t) 


application 



Figure 1: Syntax of the lambda calculus A 



A term of the form (Ax. t) is called a functional abstraction. It represents 
the function x ^ t. For example, the identity function x ^ x is written as 

(Ax. x) 

The dummy variable x here is called a bound variable, and conforms to 
the usual rules governing bound variables in mathematical formulae. For 
example, we identify expressions that differ only in the renaming of bound 
variables. 

A term of the form {t t) represents a function application. The sole means 
of computation in the lambda calculus is the operation of applying a function 
to its argument consistent with the following axiom: 

((Ax. t)v) = t [v/x] {13) 

Here v denotes a value, to be defined shortly. Reading this axiom from left 
to right defines an algebraic rewrite rule for transforming terms, substituting 
the argument v in place of the variable x into the function body. This 
transformation is called beta reduction. We will use the arrow — > to indicate 
one (and sometimes more than one) beta reduction step. A reducible term 
is called a redex. 

Unabridged lambda terms can be painful to read. For this reason, we will 
often introduce abbreviations using the symbol =. In addition, we will often 
omit parentheses according to the convention that nested lambda abstrac- 
tions associate to the right and applications associate to the left. 

Consider the simple program ((Ax.x) apple) where apple stands for some 
term in our language. With the abbreviation id = (Ax.x), this may be 
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written more legibly as (id apple), which should evaluate to apple. Indeed, 
beta reduction gives 

((Ax. a;) apple) — > apple 

in a single step. 

In general, a computation consists of a sequence of beta reductions exe- 
cuted according to some deterministic strategy until the resulting term can- 
not be reduced any further, at which point the computation terminates. 

A slightly more complicated example, which serves to show how multiple- 
argument functions can represented in terms of nested single- argument func- 
tions (a technique known as currying), is given by 

apply = A/. Xx.{f x) 
^ A/. (Ax. (/ x)) 

which represents a function that applies its first argument /, which should 
be a function, to its second argument x. Applying the identity function 
to banana should give banana. To see this, the program (apply id banana), 
which is shorthand for ((apply id) banana), is now executed by the following 
sequence of beta reductions (underlining redexes) 

((apply id) banana) = (((A/. (Ax. (/ x))) id) banana) 
— > ((Ax. (id x)) banana) 
— »• (id banana) 
— >■ banana 

Often there is more than one reducible subterm at any given step and a 
strategy is required to make the process unambiguous. For definiteness, we 
will use a call by value strategy. This works as follows: Abstractions (terms of 
the form (Ax. t)) are considered values and may not be reduced any further. 
A function application {t t) may only be reduced if both the operator and the 
operand are values. Otherwise the operator and operand must be reduced 
first. We will call the resulting calculus the call-by- value lambda calculus A^,. 
Formally, we state the syntax for values [201 1^ in figure |21 
The reduction rules are listed in figure El 
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V ::— 

X 

{\x. t) 



values: 
variable 

abstraction value 



Figure 2: Values in the call-by-value calculus A; 



ti t[ 

{h t2) {t[ t2) 

^2 ^2 

(^^1 ^2) {Vi 4) 

{\x. t) V — > t[v/x\ 



Figure 3: Reduction rules for the call-by-value calculus 

Wc will denote by {{t)) the term, when it exists, obtained by fully reducing t 
to a value. 

We will often use a less cumbersome informal notation when defining 
functions. For example, the apply function above satisfies the following prop- 
erty 

apply fx — ^{f x) 

under beta reduction. Given this specification, the translation into a lambda 
term is straightforward. 

How do we represent data in the lambda calculus? Since all we have at 
our disposal are lambda terms, we need a way of encoding data as lambda 
abstractions with specified properties. There is a technique which can be used 
for any kind of data structure, which we will illustrate with two examples: 
Natural numbers and lists. 

Let us first consider how the natural numbers may be represented. As 
with any kind of data, we need a way to construct natural numbers and a 
way to deconstruct them, extracting their constituents. 

One possible encoding is as the sequence 

0, 1 = ((succ 0)), 2 = ((succ 1)), . . . , 



{apPi) 
{app^) 
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where 



= Xx. Xy. {x id) 
succ = An. Xx. Xy. {y n) 

are the constructors.^ The above definitions were motivated by the need to 
be able to define a case expression (deconstructor) 

case of (0 — > ^2, succ m — > t^) 

which may now be taken as an abbreviation for 

^1 (A-z. ^2) (Am. ^3) 

Here z denotes a variable that does not appear free in t^. This expression 
allows us to deconstruct a natural number, extracting the ingredients that 
went into its construction, i.e., either or its predecessor m. It is indeed not 
difficult to verify the following behavior under beta reduction: 

case of (0 ^ h, succ m — > ^3) — > 
case ((succ Zq)) of (Q h, succ m — > ^3) — ^ [{{to)) / ttt] 

As an example, it is now trivial to define the predecessor function (with the 
convention that pred = 0) 

pred = Xn. case n of i - ^ - 

[ succ m ^ m 

In order to program arbitrary computations, we need to verify that the 
lambda calculus is sufficiently powerful to represent recursive functions. In- 
deed, recursion can be used to represent any kind of iterative or looping 
computation. 

^Explicitly 

= A.T. \y. (x {\w. w)) 

1 = Ax. Xy. (y Xx. Xy. {x {Xw. w))) 

2 = Ax. Ay. (y Ax. Ay. (y Ax. Ay. (x {Xw. w)))) 
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For example, in order to define addition, we need an expression add which 
behaves as follows under beta reduction: 



add m n — > case m {~ ^ 



succ k add k (succ n) 

where the subterm denoted by add on the left has copied itself into the body 
of the term on the right hand side. One of the simplest ways of achieving 
this is to define |22| 

add = (t t) (1) 

where 

t = Xf.l Am. An. case m of i " ^ r , ^^ , r A (2) 
\ I succ k ^ [J J) k (succ n) I 

In other words, t is an abstraction consisting of the body of the addition 
function with the combination (/ /) in the position where add should insert 
itself after reduction. It is a simple exercise to show that add indeed has the 
specified behavior under beta reduction. This method can be applied to any 
recursive function. 

The computation of the program (add 2 2) then proceeds via the following 
sequence of beta reductions: 

add 2 2 = add ((succ 1)) 2 

— > case ((succ 1)) of (0 2, succ k add k (succ 2)) 

— > add 1 ((succ 2)) = add ((succ 0)) 3 

— > case ((succ 0)) of (0 3, succ k add k (succ 3)) 

— > add ((succ 3)) = add 4 

— > case of (0 — 4, succ k add k (succ 4)) 

— ^4 

The above technique can be generalized to arbitrary data structures. For 
example, lists can be represented by the following constructors, which are 
entirely analogous to those of the natural numbers 

= Aa;. Xy. {x id) 
cons = Xh. At. Ax. Xy. {{y h) t) 
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where () denotes the empty hst and cons constructs a hst consisting of a 
first (head) element h followed by a list t (the tail) containing the rest of 
the elements. Again, the above definitions were motivated by the need to be 
able to define a deconstructor 

case ti of (0 t2, cons ht ^ t^) 

which may now be taken as an abbreviation for 

ti {\z. h) (A/i. At. h) 

We will often abbreviate h : t = (cons h t). We can define tuples in terms 
of lists as {xi, . . . , Xn) = Xi : X2 : ■ ■ ■ : Xn : {). Under beta reduction, we have 
the behavior 

case of (0 ^ t2, h -.t ^ t^) — > t2 
case ((<o : k)) of (() t2, h : t ^ h) ^ h [{{to)) IK {{ti))/t] 

showing how the case expression may be used to deconstruct the list, ex- 
tracting its head and tail. To see how these abstractions are used, consider 
the following recursive function 



map / list — > case list of 



0-0 

h:t^{fh): (map / t) 



which takes as input a function and a list and applies the function to each 
element of the list. The reader may verify that, for example, 

map double (4, 7, 2) — )• (8, 14, 4), double = Xx. (add x x) 

Finally, we introduce some convenient notation. Since we can represent tuples 
as lists, we can define functions on tuples using notation such as 

r - 

X{x,y). t ^ Xu.c^se u I _ ^ ^^^^ ^, f () () 

I yy -.t ^t 

For example, {X{x^y). (add x y)) (7, 7) evaluates to 14. It is also useful to 
have a notation for representing intermediate results. The let notation 

let {xi,...,Xn) = {ti,...,tn) in t 

= {X{xi, ...,Xn).t) {ti, ...,tn) 
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allows us to write terms such as 

let a; = 1 in 
let (y,2;) = (2,3) in 
add X (add y z) 

which evaluates to 6. 

3 A quantum computational model 

In this section we will construct a computational model, based on the lambda 
calculus, suitable for describing quantum computations. The language used 
will be an adaptation of the classical lambda calculus, extended with a set of 
quantum primitives. We will denote it by where the subscript stands for 
intermediate. For reasons to be discussed in the next section, this language 
is not suitable as a formal system. In particular, reduction in Aj does not 
correspond to a simple system for equational reasoning. In sectional we will 
correct these deficiencies to obtain the full quantum lambda calculus A^. 

In the classical lambda calculus, beta reduction consumes the program 
to give the result. At each step, information is discarded, which makes the 
process irreversible. For quantum computing, we need reduction rules that 
take computational states to superpositions of states in a way that is unitary 
and reversible. 

Bennett [21] showed that any classical computation can be transformed 
into a reversible computation. The construction, adapted to our situation, 
is as follows: Let x denote the term being computed, and let /? : x i— >• 
denote a single beta reduction step. Instead of the non-invertible function P, 
one considers the function x ^— {x,j3{x)), which is invertible on its range. In 
its simplest version, the computation proceeds as 

X ^ (x, (3{x)) ^ {x, (3{x),l3\x)) ^ {x, I3{x),(3\x),(3\x)) 

More complicated schemes exist that reversibly erase the intermediate steps, 
saving space at the expense of running time. Although this process does 
not end by itself, we may observe it and regard the computation as having 
terminated when (3^~^^{x) = f3^{x), at which time we may stop the machine 
by external intervention. 

Although this scheme can be used to reversibly implement computations 
in the classical lambda calculus, we will soon see that it does not work un- 
modified in the quantum case. 
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In order to represent computations involving qubits, we will add a few 
constant symbols as additional primitives to our language as in figure 01 



t ::= terms: 

X variable 

(Ax. t) abstraction 

{t t) application 

c constant 

c ::= constants: 

I 1 I I S I i?3 I cnot \ X \ Y \ Z \ ... 



Figure 4: Syntax of the intermediate language Aj 

The symbols and 1 here are primitives and should not be confused 
with the abbreviations and 1 of the previous section. Additional constants 
H, S, . . . , will denote elementary gate operations on qubits. These should 
include symbols for a universal set of elementary quantum gates jHl 1211 • 
For example, the set consisting of the Hadamard gate H, the phase gate 
S, the 71 / 8 gate R3, and the controlled- not gate cnot is universal [2^1 IHj- 
Additional primitives, such as the the Pauli gates X, Y and Z, may be 
added for convenience. 

We now allow the state of a computation to be a quantum superposition of 
terms in this language. As a model, one may imagine lambda terms encoded 
as strings of symbols on the tape of a quantum Turing machine. 

As a first example, consider an initial state written in ket notation as 

\{H 0)). 

We would like to choose the transition rules of the quantum computer in 
such a way that this string will evaluate to the Hadamard operator applied 
to |0), which should give the superposition (|0) + |1)) of the states |0) and 
|1) containing unit-length strings. The candidate reduction rule 

l(^ 0)) — -^(|0) + |1)) 

m 1)) — 4(io)-ii)) 
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is not reversible. To make it reversible, we first try the same trick as in the 
classical case 

l(^ 0)) 0);0) + |(if 0); 1)) 

^\{H 0))®i=(|0) + |l)) 

where we have factored out the common substring. The semicolon denotes 
string concatenation. In this simple example, the answer indeed factors out 
on the right. However, notice what happens if we apply this method to the 
term 

\{H {H 0))) {H 0)); {H 0)) + \{H {H 0)); {H 1))) 

-^1\{H {H 0)))® 

' {H 0); 0) + \{H 0); 1) + \{H 1); 0) - \{H 1); 1) 

Here the answer does not factor out. The fully reduced rightmost term is 
entangled with the intermediate term in the history. 

Note, however, that this scheme keeps more information than necessary. 
For reversibility, it is sufficient to record at each step only which subterm 
has been reduced, and the operation that has been applied to it. We may 
encode this in our example as follows (to be formalized below): 

\{H {H 0))) (|(_ [H _))■ {H 0)) + |(_ {H _))■ {H 1))) 



{H 0) + \{H 1) + \{H 0) - \{H 1] 

= (//_));(// J) ®|0) 

Here we have at each step replaced subterms that do not need to be recorded 
by the constant placeholder symbol _. Now the answer does indeed factor out 
on the right as required, consistent with H'^ \ 0) — |0). It is also clear that at 
each step we have kept enough information to reconstruct the previous step, 
thus ensuring reversibility. 
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V ::= values: 

X variable 

c constant 

{\x. t) abstraction value 



Figure 5: Values in the intermediate language Aj 



The computational model may now be formalized with the following rules: 
First, we extend our definition of values to include constants as in figure El 
The computational state is taken to be a quantum superposition of sequences 
of the form 

hi, . . . , hfi, t 

where hi, . . .; hn will be called the history track and t will be called the 
computational register. The classical subset of the transition rules is shown 
in figure ini 



ti — » hi, t[ 

n; (ti ta) H; (h _); (t'l is) 



^2 ^ t'2 



(apPi) 
(app^) 



H; {vi t2)-^n; (_/i2); {vi t'^) 

H; ((Ax. t) v) — > H; {{\x.tx) _); t [v/x] (Pi) if x appears free in t 

Ti; ((Ax. t) v) — > H; ((Ax._) v); t (/^s) if x not free in t 

TC; t — > Ti; _; t (Id) otherwise 



Figure 6: Operational model for the classical subset of A^ 

In these rules H denotes the (possibly empty) history track, and tx is 
obtained from t by recursively replacing all subterms that do not contain x 
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with the placeholder symbol and keeping x. More formally 



if X not free in t 



{Xy. 0, = (_. t^) 



(3) 



These rules are sufficient to make classical computations reversible, provided 
that lambda terms which differ only by renaming of bound variables have 
been identified. In this regard, we note here that in a quantum Turing 
machine model, it is possible to represent terms on the tape of the quantum 
Turing machine in an unambiguous way (e.g., using De Bruijn indices instead 
of bound variables) [HEj- 
Here is an example computation: 



|((apply id) banana)) = (( 



A/. {Xx. (/ x))) {Xz. z)) banana) 



^l(((A/. (_. 


(/ 


-))) J J 


((Ax. 


((A^.^) 


x)) banana)) 


— I(((A/. (_. 


(/ 


-))) J J 


((Ax. 


(- x)) _ 


); {{Xz. z) banana)) 


^I(((A/. (_. 


(/ 


-))) J J 


((Ax. 


L ^)) - 


); {{Xz. z) _); banana) 


— I(((A/. (_. 


(/ 


-))) J J 


((Ax. 


(- ^)) - 


); {{Xz. z) _); _; banana) 


— I(((A/. (_. 


(/ 


-))) J J 


((Ax. 


(- x)) _ 


); {{Xz.z) _); _; _; banana) 



At each step just enough information is kept to reconstruct the previous step. 
Although in this particular example, termination can be tested by observing 
and comparing the last expression in the history with _, in general we do 
not have a well-defined criterion for termination in the calculus Aj, due to 
the fact that the state may involve a superposition of several computational 
histories, some of which have terminated and others not. Thus, to observe 
termination would potentially disturb the state. This problem will be solved 
in the quantum calculus A, of sectional 

In addition, we have some extra reduction rules involving the quantum 
gate symbols such as: 

The rules for quantum primitives are summarized in figure |H1 
Here cu denotes any one of the quantum primitive symbols and U the cor- 
responding unitary transformation, while stands for or 1 in the case 
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n; {H 0))^\n; (//_))®-^(|0) + |l)) 
H; (H l))^\n; (// J) (|0) - |1)) 

Figure 7: Operational model for H 



n;{cu<P))^\n;{cuJ)(^U\<f>) {U) 



Figure 8: Operational model for the quantum primitives of Aj 

of single-bit operators, or one of (0,0), (0, 1), (1,0) or (1, 1) in the case of 
two-bit operators. For example 

\{cnot (1,0))) ^Kcnot J; (1,1)). 

4 Towards an equational theory 

While the language Aj constructed in the previous section can be used to 
describe quantum computations, reduction in Aj does not correspond to a 
simple system for equational reasoning. This makes Aj unsuitable as a formal 
proof system for quantum computation. We will discuss the problem in 
this section and resolve it in the next with the introduction of the quantum 
lambda calculus Xg. 

In the classical lambda calculus, program evaluation through beta reduc- 
tion can be regarded as a directed form of equational reasoning consistent 
with the axiom 

(Xx.t) V ^t[v/x] (P) 

Indeed, the classical lambda calculus provides both a model of computation 
and a formal system for reasoning about functions, a property we would like 
to keep in the quantum case. 

To understand the difficulty, notice what happens when a function appli- 
cation discards its argument (in other words, the argument does not appear 
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in the function body). For example 

I ((Ax. apple) banana)) — > |((Ax. _) banana); apple) 

We see that in order to maintain reversibihty, a record of the argument 
banana is kept in the history. Restricting our attention to the computational 
register, we see that its evolution is consistent with replacing the original 
expression with an equal expression according to the axiom (/5). In other 
words, in this example reduction is consistent with equational reasoning. 

However, we run into problems when the discarded subterm is in a quan- 
tum superposision with respect to the computational basis. For example, 
consider the reduction of 

\{\x. 0) {H 0)) ^ |(_ {H _)) ® i= (^\{\x. 0) 0) + \{\x. 0) 1) 



(- (H _)) ® -^(|(Ax._) 0) + I (Ax. J 1)^ ® |0) 



In the second step, a discarded subterm in a superposition is saved in the 
history and the computational register becomes |0). However, if we were to 
apply the axiom (/3) to the contents of the computational register, we would 
get the equation 



-^(^KAx.O) 0) + |(Aa;.0) 1)) = |0) 



V2' 

which is invalid since the right hand side is not a legal normalized state. 

As a second example, consider the following computation, where the inner 
function discards its argument x: 

\{{Xy.{{Xx.y) y)) {H 0))) 

^ l(- {H _))) ® (\{{Xy. {{Xx.y) y)) 0)) + \{{Xy. {{Xx.y) y)) 1)) 
1 

^71 



1 

7i 



{H _))-{{Xy.{{_.y) y))_)) 
[((Ax.O) 0)) + |((Ax.l) 1)) 
(_ {H _))-{{Xy.{{_.y) y))_))® 
\{{Xx._) 0);0) + |((Aa;._) 1); 1) 
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Now the computational register is entangled with the last expression in the 
history. Ignoring the history, the computational register would be in a mixed 
state with density matrix 

\ 
1 

However, an attempt to apply the equational axiom to the contents of 
the computational register would give 

^ (l((Aa:. 0) 0)) + |((Aa:. 1) 1))) = (1°) + l^)) 
which is clearly inconsistent. 



5 A quantum lambda calculus 

We will resolve the shortcomings of the language Aj by developing a quantum 
lambda calculus \q which has a consistent equational theory. This section 
will be somewhat heavier on the formalities, and the reader who wishes to 
see some concrete examples may wish to skip ahead to section [7| after reading 
the introductory paragraphs. 

The previous discussion suggests that the problems with equational rea- 
soning in the presence of quantum operations can be avoided by preventing 
functions from discarding arguments which may be in a superposition with 
respect to the computational basis. 

Let us call a subexpression definite with respect to the computational 
basis if it is textually the same in all branches of the superposition. For 
example, in the state 

i=(^|(Ax.O) 0) + |(Ax.O) 1) 

the subexpression (Ax. 0) is definite, whereas the argument | (|0) + |1)) is 
non-definite. Definite subexpressions may be thought of as a classical re- 
source. They can be observed without affecting the state of the computation. 
On the other hand, non-definite subexpressions represent purely quantum re- 
sources. 

To avoid the problems pointed out in the previous section, we seek a 
calculus that will keep track of whether an argument is definite or non- 
definite, and which will make it impossible to write a function that discards 
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a non-definite resource. Calculi that are resource sensitive, known as linear 
lambda calculi, have been studied intensively in recent years [221 1211 I2H1 
1^ . So-called typed linear lambda calculi are very closely related to the 
field of linear logic [201 1^- Linear logic is a resource sensitive logic where, 
for example, certain assumptions may only be used once in the course of a 
derivation. 

For our purposes it will be sufficient to study a simple untyped linear 
calculus. The syntax is a fragment of the one introduced in |2Z|, extended 
with quantum operations as in figure (HI 





terms: 


X 


variable 


{Xx. t) 


abstraction 


(tt) 


application 


c 


constant 


\t 


nonlinear term 


{Wx.t) 


nonlinear abstraction 




constants: 



I 1 I I 5 I i?3 I cnot \X\Y \Z 



Figure 9: Syntax of the quantum calculus \q 

Here terms of the form \t are called nonlinear. Nonlinear terms will be 
guaranteed to be definite with respect to the computational basis, and may 
be thought of physically as classical strings of symbols that may be discarded 
and duplicated at will. On the other hand, linear terms may be non-definite, 
possibly containing embedded qubits in superpositions with respect to the 
computational basis. Abstractions of the form (A!x. t) denote functions of 
nonlinear arguments. In an abstraction of the form (Ax. i), the argument is 
called linear. 

A functional abstraction may use a nonlinear argument any number of 
times in its body, or not at all. On the other hand, a linear argument must 
appear exactly once in the function body (hence the name linear). 

To enforce these rules, we require terms to be well-formed. This corre- 
sponds to the constraint that linear arguments appear linearly in a function 
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body, and that all free variables appearing in a term \t refer to nonlinear 
variables j2Hl- In the following examples, the terms in the left column are 
well-formed, while those in the right column are ill-formed ^ 



Well-formedness is a property that can be checked syntactically. For com- 
pleteness, we formally state the rules for well-formedness |28j, which the 
reader satisfied with the above informal characterization may skip, in figure 



These rules may be related to the typed linear calculi described in PHIITF] 
by erasing the type annotations from the typing rules of the latter. Here F 
and A denote contexts, which are sets containing linearity assumptions of 
the form x and \x, where each variable x is distinct. If F and A are contexts 
with no variables in common, then F, A denotes their union. For example, 
the rule —o-E implicitly assumes that F fl A = 0. Rules may be read as 
follows: For example, the promotion rule says that if t is a well-formed term 
under the assumption that xi to nonlinear, then \t is a well-formed 

term under the same assumption. The condition Ffl A = in (^-E) ensures 
that a linear variable can only appear once in the body of a formula. The 
weakening and (^-I) rules allow a function to discard a nonlinear argument, 
whereas the contraction and (— o-E) rules allow us to duplicate a nonlinear 
argument any number of times in the body of a function. 

The well-formedness constraint prevents us from writing a function which 
discards a linear argument. However, this is not sufficient to prevent unsafe 
computations without further specification of the substitution order. To see 
this, consider the expression ((A!x. 0) \{H 0)), which is well-formed. The 
problem is that we are allowed to use ! to promote the expression {H 0) to a 
nonlinear value, which can then be discarded. If we were allowed to reduce 

^Notice that while weh-formedness guarantees that hnear resources will be used ap- 
propriately, it does not guarantee that terms are meaningful. For example the term 
{Xy. (A !z.O) y) is well-formed, but may or may not get stuck at run-time, according to 
the operational model of figure El when applied to a linear or nonlinear argument respec- 
tively. A typed calculus would be needed to specify which terms can be legally substituted 
for y. For recent progress in this direction, see |16| . 



(A!x.O) 



(Ax. 0) 





{X\x. {x x)) 

{Xy. (Ab.y)) 

{X\y.\{X\x.y)) 



(Ax. (x x)) 
{Xy. {Xx.y)) 
{Xy.\{X\x.y)) 
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he 

x\- X 
\x\. . . . , \x^ h t 
. . . , ^'Xfi \~\t 

r,x h t 

r,!a: h t 

r, \x, \y h t 
r, l^; h t [2;/a;,2;/|/] 

r ht 

r,!x h t 

r,a; h t 

rh (Ax.t) 

r,!a: ht 

rh (Ab.i) 

r h A h i2 

r, A h (ti t2) 

Figure 10: Rules for well-formed terms in the quantum calculus 

the subterm [H 0) first, equational reasoning would give 

|((A!x.O) \{H 0))) = i=(^|((A!x.O) !0)) + |((A!x.O) !l))) 
= ^/2|0) 

which is an invalid equation since the last line is not a valid normalized state. 
On the other hand, if we consider \{H 0) as an irreducible value, we may use 
beta reduction immediately to obtain 

|((Ab.O) \{H 0))) = |0) 

which is a valid result, since we are discarding the unevaluated expression 
\{H 0), which is definite. 



Const 
Id 

Pomotion 
Dereliction 
Contraction 
Weakening 
^-I 
^-I 
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To prevent terms of the form \ t from being evaluated, we follow Abramsky 
j2ni and extend our definition of values as in figure [TT] 



V :: = 


values: 


X 


variable 


c 


constant 


(Ax. t) 


linear abstraction 


(A!x. t) 


nonlinear abstraction 


H 


[-suspension 



Figure 11: Values in the quantum calculus Xq 



The computational model is described in figure [121'^, where t is defined as in 



ti — » hi, t'l , . 

n; (ti t2)^n; (hi J; {t[ t^) ^™ 



t2 > h2] t'2 



(«m) 



{vi t2) H; (_ /is); (^^1 t'2) 

H; ((Ax. t) v) H; ((Ax. t,) _); t [v/x] (/3) 

H; ((A!x. t) H') — > H; ((Alx.I^;) _); t [t'/x] (!/3i) if x appears free in t 

H; ((A!x. t) W) — > H; ((A!x._) W); t (!/?2) if x not free in t 

\n; {cu <p)) \n; {cu_))^U\<P) (U) 

H; t — > H; _; t (Id) otherwise 



Figure 12: Operational model for the quantum lambda calculus Xq 

According to these rules, quantum superpositions can only be created by 
evaluating terms containing quantum primitives. The result of applying a 

•^See on] for related operational interpretations of linear lambda calculi. Our 

evaluation model recomputes l-closures (see |H2p. 
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quantum gate is a linear value, not preceded by a !. As we prove below, 
there is no way of including such a linear value in a nonlinear subterm. 
It follows that subterms that may be quantum non-definite will never be 
discarded since, by and (!/32), nonlinear functions can only be applied 
to nonlinear terms. 

Note that when a nonlinear function encounters a linear argument, it 
simply gets stuck. More precisely, the rule (Id) applies. 

The above reduction rules may create superpositions. However, such 
superpositions are not arbitrary. Indeed, terms in a superposition will only 
differ in positions containing the constants and 1. Otherwise they will have 
the same shape. We may formalize this by defining two terms to be congruent 
if they coincide symbol by symbol except possibly in positions containing 
or 1. It then follows that 

Lemma 5.1. All terms in a superposition obtained via a reduction sequence 
from a definite initial term are congruent. 

Proof. The proof is by a simple induction on the length of the reduction 
sequence, analyzing the reduction rules case by case. □ 

Another case by case induction argument may be used to prove that 
reduction preserves well-formedness. More precisely 

Lemma 5.2. Ift is well-formed and |7Y; t) — > I'^ii ^i)? ^^^'^ terms 

t[ appearing in the resulting superposition are well-formed.'^ 

Because terms appearing in a superposition have the same shape, it makes 
sense to talk about specific subterms of the expression in the computational 
register. We can therefore formulate the following lemma: 

Lemma 5.3. Starting from a definite initial term, any \-suspension subterm 
occurring during reduction is definite with respect to the computational basis. 

Proof. This follows by induction on the length of the reduction sequence. 
The initial term is definite by assumption. Assume that the lemma holds 
after n steps. We have argued that all terms in a superposition obtained 

from a definite initial term are congruent. They therefore have the same 

"^Thanks to one of the referees for suggesting improvements in the exposition of this 
section. 
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structure of subterms and the same reduction rule applies to them all. Since 
we have argued that these terms are well-formed, there are then three ways 
in which we may obtain a !-suspension subterm after n + 1 steps. First, 
the suspension may not be part of the redex, in which case it is included 
unmodified in the resulting expression. Second, it may be the result of beta 
reduction of an application of the form 

{Xx.{---x---)) {...\t---) 

where \t is definite by the induction assumption. The result is (■■■(■■■ !t ■■■)■•• ), 
where \t has been copied without modification. Third, it may be the result 
of beta reduction of an application of the form 

{X\x. ... \{...x ■■■)■■■) \t 

where It and \{- ■ - x ■ ■ ■) are definite by the induction assumption. This cre- 
ates a suspension !(• • - i • • •), which is definite because all its subterms are 
definite. This completes the proof. □ 

It is worth pointing out that we cannot create possibly non-definite sus- 
pensions by reducing terms like 

{Xx. ...\{...x---)---) {H 0) 

because x is linear, which implies that \{- ■ • x • ■ •) is not a well- formed sub- 
term. 

Lemma 5.4. Given a definite initial term, the contents of the history track 
remains definite throughout reduction. 

Proof. We have argued that all terms in a superposition obtained from a 
definite initial term are congruent. They therefore have the same structure of 
subterms and the same reduction rule applies to them all. Since our reduction 
rules allow only !-suspensions to be discarded (which saves a copy in the 
history), and since !-suspensions are always definite by the previous lemma, 
the result follows by induction on the length of the reduction sequence. □ 

Since both the history and the shape of the term in the computational 
register remain definite throughout reduction, we can state the following 
conclusion 
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Corollary 5.5. Termination can be tested without disturbing the computa- 
tion by observing the last term in the history. When this term becomes equal 
to the placeholder _, the result can be read off from the computational register. 

Not only is termination a "classical" property, but so is the entire shape 
of the term itself, i.e., term shapes can be implemented in classical memory, 
and only their data slots need to point to qubits on a quantum device.^ 

The fact that the history remains definite in Xg eliminates the specific 
impediments to setting up an equational theory that were pointed out in the 
previous section. Indeed, since the state of the computation is now always 
guaranteed to be a direct product |7i) ® |c) of the history \TC) and the com- 
putational register |c), reduction can never lead to a computational register 
|c) that is in a mixed state. In addition, since remains definite, the 
restriction of the reduction rules to the computational register will preserve 
the normalization. We are therefore led to the following theorem: 



h - 




(apPi) 


ih t2) - 


- it'l t2) 


t2 - 




(app^) 


K t2) - 


- (^1 t',) 


(Ax. t) V — 


t [v/x] 




(A!x. t) \t' - 


t [t'/x] 


m 


\CU 0) - 


-f/|0) 





Figure 13: Reduction rules for the quantum calculus Xg 

Theorem 5.6. In the quantum calculus Xg, the evolution of the computa- 
tional register is governed by the reduction rules of fiaureM'Ji 

Proof. This easily follows from a case-by-case analysis of the computational 
rules of figure [T21 

would like to thank one of the referees for suggesting this improved formulation. 
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For example, consider the rule (!/32) applied to a state of the form \H) ®\c) 
where |c) is a normalized superposition of the form 

Y,C^ \{■■■^ {{Wx.k) \t') ■■■,)) 

i 

in which, by lemma IK?T| all terms have the same structure and by lemma 
the subterm \t' does not depend on i. This then reduces to 

i 

= \n; (_((A!x._) \t')_)) ® |(- ■ [t'/x] ■ ■ -,)) 

i 

with the computational register in the normalized state q |(- ■ -j ti [t'/x] ■ ■ -i)), 
consistent with applying the reduction rule {X\x. ti) \t' — > ti [t'/x] to its con- 
tents. □ 

The big win is that we now have a simple set of reduction rules that can 
be used to reason about the computation without having to keep track of 
the history. 

In order to define an equational theory for this calculus, we will simply 
define a notion of equality that is compatible with the the reduction rules of 
figure Uni Intuitively, reduction should be understood as a simple algebraic 
operation of replacing subterms with equal subterms. However, we need to 
take into account that reduction may not take place inside ! -suspensions. 

We therefore need to introduce algebraic rules governing just when we 
can replace subterms in an expression with equal subterms [^|^. One way 
to do that is to introduce the notion of a term context, which are expressions 
with a hole [ ] in place of a subexpression: 

C,Ci::=[]\{tC)\ (C t) I (Ax. C) I (A!a;. C) \ ^ q d, 

i 

where the last term denotes a superposition of shape-congruent contexts. It 
is important to note that there are no contexts of the form !(■■■[]■■■). As 
a result, subterms preceded by ! will be opaque in the sense that we will not 
be able to perform substitutions under the ! sign. 
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Definition 5.7. The equational theory of is the least equivalence relation 
= containing the reduction relation ( — >) of figure El and which is closed 
under substitution in term contexts fl^ |^ . In other words 

{subst) 



C[h] = C[h] 



where C is an arbitrary term context and C[t] denotes the textual replace- 
ment of the hole in C by the term t, extended by linearity to superpositions 
of congruent terms and congruent contexts, i.e., for t = ^^Citi, we define 
C[t] = Ei Ci C[U] and for C = q we define C[t] = ^. q a[t]. 

An alternative way of presenting the equational theory is by listing a set 
of axioms and rules of inference as in figure IT^ These rules should again be 



t = t {refl) 

ti = h 



t2 = tl 
ti = t2 t2 = ts 



tl ■■ 


= t3 


tl=t2 


ts = t^ 


{tl h) 


= {t2 ti) 


tl 


= t2 


Xx. tl 


= Xx. t2 


tl 


= t2 


Xlx.ti 


= X\x.t2 



(sym) 
(trans) 
(app) 

(Ai) 
(A2) 



{Xx. t) v = t [v/x] {(3) 
{X\x. t) \t' = t [t'/x] {\p) 

|(cc;0)) = f/|0) {U) 



Figure 14: Equational proof system for the quantum calculus Xq 

understood as extending via linearity to congruent superpositions of terms. 
In this formulation, the rules {app), (Ai) and (A2) are together equivalent to 
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the term context substitution rule (subst) above. Again, there is no rule that 
permits substitutions inside !-suspensions. 

Theorem 5.8. In the quantum lambda calculus, the evolution of the compu- 
tational register proceeds by replacing terms by equal terms according to the 
equational theory of definition \5 . 7[ 

Proof. True by construction. □ 



6 Recursion and a fixed point operator 

Recursive functions may be defined in the calculus Ag in a way analogous 
to that described in section |21 We simply replace {t t) in equation ((T)) with 
{t where now t = A!/. (■■■(/!/)■•■)• 

Here we describe a related approach based on so-called fixed point com- 
binators. A fixed point operator suitable for the linear lambda calculus is 
given by the following adaptation of the classical Turing combinator 

fix ^ ( (A!«. A!/. (/!((« \u) !/))) 
!(A!«.A!/. (/!((« \u) !/)))) 

It is easy to check that under reduction 

fix \t — >t !(fix \t) 

where the !-suspension prevents further reduction of the term in brackets. 
Recursive functions can be defined as follows: If 

t = Wf.u 

then it easily follows that 

fix \t — >u[{\\x \t)/f] 

In other words, fix \t copies itself into the body u oi t under reduction, as 
required for recursion. 
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7 Examples of algorithms 



We are now ready to formulate some algorithms in the quantum lambda 
calculus. First, we reproduce some classical constructions, now decorated 
with the proper nonlinearity annotations. 

First, we introduce list constructors that will enable us to build lists of 
linear values (qubits or structures containing qubits) 

= Wx.Wy. {x id) 
cons = \h. \t. \\x. \\y. {(y h) t) 

with abbreviations h : t = (cons h t) and (xi, . . . , x„) = Xi : X2 Xn '■ 

as before. Since the arguments \x and \y above are nonlinear, we need to 
redefine our case abbreviation as follows: 



case ^1 of (0 — > ^2, ^ : ^ — ^ h) 



now stands for 



h \{X\z.h) \{\h.Xt.h) 
Deutsch's algorithm E] can be very simply expressed as follows: 

deutsch Uf — > let {x,y) = Uf {{H 0), {H 1)) in 

{{H x),y) 

Here the argument Uf is assumed to be a function that takes (x, y) to 



|0) 
|1) 



H 



H 



Uf 



H 



Figure 15: Deutsch's algorithm 

(x, ?/ © f{x)), where / is some (unknown) function of one bit. For example, 
if / is the identity function, then we should take Uf to be cnot . Indeed, the 
reader may check that 



|deutsch cnot) 



|i)®^(|o)-|i) 
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where the first bit 1 = /(O) © /(I) indicates that the function is balanced, 
as required. 

Let us write a simple expression that creates an EPR pair 

epr = cnot {{H 0),0) 

The quantum teleportation gate array with deferred measurement |3] can 
easily be translated into the following code: We create an EPR pair and 
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Figure 16: Quantum teleportation 

pass the first EPR qubit, along with the unknown qubit x to be teleported, 
to Alice. The outcome {x' ,y') of Alice's computation then gets sent to Bob, 
who has access to the second EPR qubit 62 

teleport x — > let (ei, 62) = epr in 

let {x',y') = alice (x, ei) in 
bob {x',y',e2) 

Here 

alice {x,ei) — ^ let {x',y') = cnot (x,ei) in {{H x'),y') 

and 

bob {x',y\e2) — > let {y",e'2) = cX (^',62) in 
let (x", e'a) = cZ {x', e'^) in 
[x ,y ,62) 

The outcome of the computation consists of the list of three qubits (x", y", e'g)- 
The teleported qubit is 63, but notice how linearity requires us to keep the 
other two qubits in the answer. The reader may check that throughout 
the computation, linear arguments are used exactly once. Implementing the 
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conditional operations cX and cZ in terms of the primitive constants is left 

as an easy exercise. 

Given recursion and lists, the map function, which applies a given function 
/ to each element of a list, may be defined as 



map !/ list — case list of 



0-0 

h:t^{fh): (map !/ t) 



The arguments list, h and t may refer to qubits or data structures containing 
qubits and are therefore chosen linear. The expression is well-formed because 
list, h and t are each used exactly once. 

It is now trivial to define a program that computes a uniform superposi- 
tion of a list of qubits by applying the Hadamard gate to each qubit in the 
list: 

list — ^ map \H list 

For example, we may evaluate 

\H^- (0, 0)> \ (^1(0, 0)) + 1(0, 1)) + 1(1, 0)) + 1(1, 1)) 

Note that the well-formedness conditions may be somewhat subtle, as the 
following example illustrates. A naive attempt at defining an append function 
that concatenates two linear lists 



append x y — > case x of 



0-2/ 

h:t^h: (append t y) 



fails to be well-formed. The problem can be seen by expanding the case 
abbreviation 

X \{\z. y) \{Xh. Xt. {h : (append t y))) 

Since y is a linear variable, we may not promote the A subterms to nonlinear 
values with the prepended !. An alternative definition that does work is 

^ {Xu. u) 

h : t ^ Xu. {h : (append t u)) 



append x y — > ^case x of 

Next we define a reverse function 
reverse list — > case list of 



0-0 

h : t ^ append (reverse t) (h) 
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Figure 17: The quantum Fourier transform (without reversal) 

The quantum Fourier transform [HHl ESI ISZj can now be defined as a direct 
translation of the corresponding quantum circuit [S] as follows: 

fourier list — * reverse fourier' list 

where 

r - 

fourier' list — > case list of < h : t ^ let h' : t' = phases {H h) t !2 in 

( h' : (fourier' t') 

recursively applies the appropriate conditional phase operations to the first 
qubit in the list, using the helper function 

phases target controls \n 

' ^ (target) 



case controls of < 



control: t — >■ let [control', target') 

= [cR \n) [control, target) in 
let target" : t' 

= phases target' t !(succ n) in 
target" : control' : t' 

Here [cR \n) composes an appropriate combination of elementary gates to 
implement a conditional phase operation with phase 27ri/2"'. Since this is 
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essentially a classical computation and depends on the particular set of prim- 
itive constants chosen, we will not write it out here. 

Note that we have assumed that the classical construction of the natural 
numbers may be adapted to the quantum lambda calculus. That this is 
possible for all classical constructions follows from the fact that there is an 
embedding of the classical lambda calculus into the linear lambda calculus, 
as shown in formula (jU) of section |H1 

8 Relating Ag to quantum Turing machines 

In this section we will sketch a proof of the following theorem, leaving a more 
rigorous analysis to future work: 

Theorem 8.1. The computational model provided by the quantum lambda 
calculus \q is equivalent to the quantum Turing machine. 

Proof. First, we argue that the quantum lambda calculus Ag may be effi- 
ciently simulated on a quantum Turing machine. 

In \q the current state of the computation consists of a superposition 
of term sequences of the form ?i; t, which may be encoded as strings of 
symbols on the tape of the quantum Turing machine. By lemma I^TTl term 
sequences in different branches of the superposition are congruent, and the 
same reduction rule will apply for all branches at each time step. The subset 
of \q not involving quantum operations consists of a set of reversible classical 
rewritings, which can be unitarily and efficiently implemented on a quantum 
Turing machine by |21 EH EHj ■ The fragment involving quantum operators 
again involves simple classical rewritings followed by a unitary transformation 
involving one or two symbols on the tape. Once again, the methods of 
[21 EHl EH] may be used to construct a quantum Turing machine that can 
execute these transformations. This completes the proof of the first half of 
the equivalence. 

Next we argue that a quantum Turing machine can be efficiently simu- 
lated by the calculus Ag. 

Yao shows in |^ that for any quantum Turing machine T, there is quan- 
tum circuit Cn,t that efficiently simulates T on inputs of size n after t steps. 
The circuit family Cn,t may be efficiently constructed via a classical compu- 
tation. But Ag is universal for classical computation. This follows from the 
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fact that the classical call-by-value lambda calculus may be embedded in 
via the following translation, adapted from j2H] 



(ti t2r = {{x\z.z) ti) t; 

X* = \x (4) 

{Xx.ty = \{x\x.t*) 

So, given the specification of a quantum Turing machine and an input of 
length n, a classical computation in Ag first constructs a representation of 
the appropriate quantum circuit family C„ It then follows the circuit di- 
agram and applies the appropriate quantum operations one by one to the 
input. Since Xq has primitive quantum operations available corresponding to 
a universal set of quantum gates, this proves the second half of the equiva- 
lence. □ 



9 Related work 

In a series of papers, Henry Baker [101 UTJ H^I develops an untyped linear 
language based on Lisp. His language is similar to the classical fragment of 
the lambda calculus developed in the current article. It served as the initial 
inspiration for the linear approach followed here. 

Ideas stemming from linear logic have been used previously by Abramsky 
in the study of classical reversible computation [IHj . 

One of the earlier attempts at formulating a language for quantum com- 
putation was Greg Baker's Qgol jH]. Its implementation (which remained 
incomplete) used so-called uniqueness types (similar but not identical to our 
linear variables) for quantum objects The language is not universal for 
quantum computation. 

The language QCL, developed by Omer, is described in [IHl El- QCL 
is an imperative language with classical control structures combined with 
special operations on quantum registers. It provides facilities for inverting 
quantum functions and for scratch space management. No formal program 
calculus is provided. A simulator is publicly available. 

Another imperative language, based on C++, is the Q language developed 
by Bettelli, Calarco and Serafini 48 . As in the case of QCL, no formal 
calculus is provided. A simulator is also available. 
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A more theoretical approach is taken by Sehnger in his description of the 
functional language QPL This language has both a graphical and a 

textual representation. A formal semantics is provided. 

The imperative language qGCL, developed by Sanders and Zuliani [^O] . 
is based on Dijkstra's guarded command language. It has a formal semantics 
and proof system. 

A previous attempt to construct a lambda calculus for quantum com- 
putation is described by Maymin in [5T]. However, his calculus appears to 
be strictly stronger than the quantum Turing machine [22] • It seems to go 
beyond quantum mechanics in that it does not appear to have a unitary 
and reversible operational model, instead relying on a more general class of 
transformations. It is an open question whether the calculus is physically 
realizable. 

A seminar by Wehr suggests that linear logic may be useful in con- 
structing a calculus for quantum computation within the mathematical frame- 
work of Chu spaces. However, the author stops short of developing such a 
calculus. 

Abramsky and Coecke describe a realization of a model of multiplicative 
linear logic via the quantum processes of entangling and de-entangling by 
means of typed projectors. They briefly discuss how these processes can be 
represented as terms of an affine lambda calculus i 54j . 

10 Conclusion 

In this article we developed a lambda calculus \q suitable for expressing and 
reasoning about quantum algorithms. We discussed both its computational 
model and its equational proof system. We argued that the resulting calculus 
provides a computational model equivalent to the quantum Turing machine 
and is therefore universal for quantum computation. 

There are many possible directions for future work. The proof of Turing 
equivalence should be fleshed out. Formal issues relating to consistency and 
semantics need to be addressed further. While our computational model 
provides an operational semantics, the problem of providing a denotational 
semantics is open. The formalism of jHS] may be useful in this regard. 

In this article, the introduction of a linear calculus was motivated by re- 
quiring consistency of its operational model with equational reasoning. The 
fact that linear arguments, denoting quantum resources, may not be dupli- 
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cated suggests a separate motivation for linearity, not addressed here, based 
on the no-cloning theorem |^ . 

While our calculus is untyped, it would be interesting to investigate typed 
linear calculi with quantum primitives and, via the Curry-Howard correspon- 
dence, the corresponding generalizations of linear logic [SHI 13^. We might 
mention that there have been prior attempts to relate linear logic to quantum 
mechanics, starting with a suggestion by Girard pUl lUUl IHT|. 

On the practical side, the calculi described in this paper may be used as 
a programming language for prototyping quantum algorithms. Indeed, the 
algorithms exhibited in this article were transcribed into Scheme for testing. 
The simulator, which was also written in Scheme, is available upon request 
from the author. 

It is our hope that the field of quantum computation, like its classical 
counterpart, may benefit from the insights provided by the alternative com- 
putational model provided by the quantum lambda calculus. 

Note added 

Since the first version of this paper was written, some progress has been made 
by the author in devising a typed version, with accompanying denotational 
semantics, of a fragment of the quantum calculus described here \16i . 
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